The effect of modified gravity on weak lensing 
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We study the effect of modified gravity on weak lensing in a class of scalar-tensor theory that 
includes f{R) gravity as a special case. These models are designed to satisfy local gravity constraints 
by having a large scalar-field mass in a region of high curvature. Matter density perturbations in 
these models are enhanced at small redshifts because of the presence of a coupling Q that charac- 
terizes the strength between dark energy and non-relativistic matter. We compute a convergence 
power spectrum of weak lensing numerically and show that the spectral index and the amplitude of 
the spectrum in the linear regime can be significantly modified compared to the ACDM model for 
large values of of the order of unity. Thus weak lensing provides a powerful tool to constrain 
such large coupling scalar-tensor models including f{R) gravity. 



I. INTRODUCTION 

The observations of the Supernovae la (SN la) in 1998 
[l| opened up a new research paradigm known as Dark 
Energy (DE). In spite of the tremendous effort over the 
past ten years, we have not yet identified the origin of 
DE responsible for the late-time accelerated expansion. 
Many DE models have been proposed so far to allevi- 
ate the theoretical problem of the cosmological constant 
scenario 0, We can broadly classify these models 
into two classes: (i) "changing gravity" models and (ii) 
"changing matter" models. The first class includes f{R) 
gravity scalar-tensor models Q and braneworld mod- 
els whereas scalar-field models such as quintessence 
01 and k-essence ^ are categorized in the second class. 

While changing matter models lead to dynamical evo- 
lution for the equation of state of DE, it is not easy to 
distinguish them from the cosmological constant scenario 
in current observations. Meanwhile, if we change gravity 
from General Relativity, the models need to pass local 
gravity tests as well as cosmological constraints. In this 
sense it is possible to place stringent experimental and 
observational constraints on changing gravity models. 

In fact there have been a burst of activities to search 
for viable modified gravity DE models. In the so-called 
f{R) gravity where / is a function of the Ricci scalar R, 
it was found that the model f{R) = R — a/R" {a, n > 0) 
proposed in Refs. Q is unable to satisfy the stability 
condition {J^rr = d^f /dR^ > 0) for perturbations 0, 
cosmological viability [lO| and local gravity constraints 
(LGC) [ll|. Recently a number of authors proposed vi- 
able /(^l DE models that satisfy all these requirements 
[13, El, Ei, [H, [H, [13, [H, [H, [lOl . For example, the model 
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/(i?) = R—aR" with a>0,0<n<lis consistent with 
LGC for n < 10~^° [2l[ while at the same time satisfying 
stability and cosmological constraints. However it is diffi- 
cult to distinguish this model from the ACDM cosmology 
because of the tight bound on the power n coming from 
LGC. 

The f{R) models proposed by Hu and Sawicki [l^ and 
Starobinsky [l^ are designed to satisfy LGC in the re- 
gion of high density where local gravity experiments are 
carried out. Moreover it is possible to find an appre- 
ciable deviation from the ACDM model as the Universe 
evolves from the matter-dominated epoch to the late- 
time accelerated era. In fact the equation of state of DE 
in these models exhibits peculiar evolution at small red- 
shifts [Tsj . In addition, for the redshift smaller than 
a critical value Zk, the growth rate of matter density per- 
turbations is larger than in the case of General Relativity 

MM- 

Recently the analysis in f{R) gravity was extended to a 
class of scalar-tensor DE models, i.e., Brans-Dicke theory 
with a scalar field potential V{cj)) [22|. By introducing a 
constant Q with the relation 1/(2(5^) = 3 -|- 2ujb'd {^^bu 
is a Brans-Dicke parameter), one can reduce this the- 
ory to the one given by the action ([2]). The constant Q 
characterizes the coupling between dark energy and non- 
relativistic matter. If the scalar field (f) is nearly massless, 
the couphng is constrained to be |Q| < 10"'^ from solar 
system experiments [l^]. However, if the field is mas- 
sive in the region of high density, it is possible to satisfy 
LGC even when \Q\ is of the order of unity. In fact, in the 
context of gravity {Q = — 1/\/6), the models of Hu 
and Sawicki 15 1 and Starobinsky [l5| are designed in such 
a way that the field is sufficiently massive in the regime 
R^ Ro (i?o is the present cosmological Ricci scalar) and 
that the mass becomes lighter as R approaches Rq. For 
general coupling Q, the potential given in Eq. ([5]) can 
be compatible with both local gravity and cosmological 
constraints. 
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The scalar-tensor models mentioned above show devi- 
ations from the ACDM model at late times and hence 
they can leave a number of interesting observational sig- 
natures. In Ref. [2^ several bounds on the coupling Q 
and model parameters were derived by considering the 
evolution of matter density perturbations as well as LGC. 
It was found that there exists allowed parameter space of 
model parameters even when \Q\ is of the order of unity. 

In this paper we shall study the effect of such modified 
gravity models on weak lensing observations [2^. Since 
weak lensing carries the information of perturbations at 
low redshifts, it is expected that this sheds light on reveal- 
ing the nature of DE [H, [H, [H, [13, [H, ^ [10, [3l|, IH] . 
In Refs. m, [13 a convergence power spectrum of weak 
lensing was derived in scalar-tensor theories with the La- 
grangian density £ = F{(j))R/2 - (V0)V2 - V{(j)). In 
these theories a deflecting lensing potential $wi is modi- 
fied compared to General Relativity due to the different 
evolution of gravitational potentials. This gives rise to 
the change of the convergence power spectrum, which 
provides a powerful tool to distinguish modified gravity 
from the ACDM model. 

The lensing potential $wi is sourced by matter density 
perturbations. The equation for matter perturbations in 
scalar-tensor models was derived in Ref. jUj under the 
approximation on sub- horizon scales (see also Ref. [36j). 
This analysis can be generalized to the theories with the 
Lagrangian density f{R,(f>,X) (where X = — (V0)^/2), 
in which <I>wi was obtained analytically [13 • The DGP 
braneworld model also leads to the modification to the 
lensing potential [11] . Thus the effect of modified gravity 
generally manifests itself in weak lensing observations. 

In this work we focus on scalar-tensor models ^ with 
a large coupling Q and evaluate the convergence power 
spectrum to find signatures of the modification of gravity 
in weak lensing. This analysis is general in the sense that 
f{R) gravity is included as a special case. In Sec. [Til we 
review our scalar-tensor models and present cosmological 
background equations to find dark energy dynamics. In 
Sec. lIIII we derive the form of the convergence power spec- 
trum as well as the equation for the deflecting potential 
^wi- In Sec. [IVj we compute the convergence spectrum 
numerically and estimate the effect of modified gravity 
on weak lensing. We conclude in Sec. [Vj 



II. MODIFIED GRAVITY MODELS 



We start with the following action 



S 



1 



(1) 



where x is a scalar field coupled to the Ricci scalar R, 
wbd is a constant parameter, V{x) is a field potential, 
and Sm is a matter action that depends on the metric 
Qf^i, and matter fields ^m- The action ([1]) corresponds 




to Brans-Dickc theory [1^ with a potential V{x)- In the 
following we use the unit 8ttG = I, but we restore the 
bare gravitational constant G when it is required. 

Setting X = F = e~^'^'^, where Q is a constant and 
(j) = — l/(2Q)lnx is a new scalar field, we find that the 
action llj is equivalent to 



^FR-^{l-6Q')F{S/<j>f -V 
), (2) 



where Q is related with the Brans-Dicke parameter wbd 
via the relation I/(2Q2) = 3 + 2wbd The f{R) 

gravity corresponds to the coupling Q = —l/\/6, i.e., 
C^BD = [13. 

In the absence of the potential V the Brans-Dicke 
parameter is constrained to be wbd > 4.0 x 10^ from 
solar system experiments (4]| . which gives the bound 
\Q\ < 2.5 X 10"'^. If the potential V is present, it is 
possible to satisfy solar system constraints even when 
IQI is of the order of unity by having a large mass in a 
high-curvature region. In the context of /(i?) gravity, 
the following model is designed to satisfy LGC [l8j : 



fiR) = R-M^-iR/Rc)-n 



(3) 



where ^, i?c, n are positive constants, and Rc is roughly 
of the order of the present cosmological Ricci scalar Rq. 
Note that this satisfies the stability condition J^rr > 
for R > Ri is a Ricci scalar at a late-time de-Sitter 
point) imlikc the model f{R) = R-a/R"- {a,n > 0) [H, 
[l8j . In the limit R^ Rc the above model approaches the 
ACDM model, which allows a possibility to be consistent 
with LGC in the region of high density. 

In fact, the model © satisfies LGC for n > 0.9 [U 
through a chameleon mechanism (4^ because of the pres- 
ence of an effecive potential V = {RF — f)/2 with the 
dynamical field = (V6/2)lni^. The field potential in 
this case is given by 



^Rc 



1 



2n + l 



(2?l^)2"/(2" 



+1 



■(1 



(4) 



The models proposed by Hu and Sawicki 15[ and by 
Starobinsky [1^ reduce to this form of the potential 
in the high-curvature region {R ^ Rc) where local 
gravity experiments are carried out. When R ^ Rc 
the field (f> is almost frozen at instantaneous minima 
around (f> ~ characterized by the condition e^'*/^ = 
1 — 27i^(i?/i?c)^^^"^^' with a large mass squared = 

2Ti-j-2 

Vtj,tj, cx (/)~ , These minima are sustained by an effec- 
tive coupling Q between non-relativistic matter and the 
field (t> U^. 

For arbitrary coupling Q with the action ([2]), one can 
also construct viable models by generalizing the analysis 
of f{R) gravity. An explicit example of the potential 
consistent with LGC is given by 0| 



V{cb) =Vi[l- C{1 



(5) 
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where Vi > 0, C > 0, < p < 1. This is motivated 
by the potential which means that the f{R) model 
Q is recovered by setting p = 2n/{2n +1). The anal- 
ysis using the potential ([5]) with the action ^ is suffi- 
ciently general to understand essential features of mod- 
ifed gravity models that satisfy local gravity and cosmo- 
logical constraints. As p gets closer to 1, the field mass in 
the region of high-curvature tends to be heavier so that 
the models are consistent with LGC. In Ref. [l^ it was 
found that the constraints coming from solar system tests 
and the violation of equivalence principle give the bounds 
p> 1-5/(9.6 -IniolQI) andp> 1-5/(13.8 -IniolQI), 
respectively. In /(i?) gravity with the potential ([4]) these 
bounds translate into the conditions n > 0.5 and n > 0.9, 
respectively pll |. 

Let us review cosmological dynamics for the action ^ 
with the potential ([5]) in the flat Friedmann-Lemaitre- 
Robertson-Walker (FLRW) metric, ds^ = -df + 
a^(t)dx^, where t is cosmic time and a{t) is the scale 
factor. As a source term for the matter action Sm, we 
take into account a non-relativistic fluid with energy den- 
sity Pm and a radiation with energy density p^d- These 
obey the usual conservation equations pm + 3Hp„i ~ 
and prad + 4_ff/9rad = 0, whcrc H = d/a. The variation of 
the action ([2]) leads to the following equations of motion: 

3FH^ = i(l - 6Q2)F02 ^jjp ^ ^ p^^^ ^ 



The effective equation of the system is defined by 

U-eff = -1 - 2ij/(3i/2) . (15) 

Using Eqs. (I6|)-(l8|), we obtain the following equations 

d.Ti 
A/V 



— {\X2 - V6xi) + ^— 



(5-6Q2)a;2 



-2\/QQxi — 3x2 



da:2 
d/V 

dX3 

dN 
where N 



f(2Q 



1 



H 



XI j^, (16) 



(17) 
(18) 



In (a) is the number of e-foldings. We note 
that the variable F satisfies the equation of motion: 
dF/dN = ^2V6QxiF. 

There exists a radiation fixed point: (xi, X2, x^) = 
(0,0,1) for this system. During radiation and matter 
eras, the field (j) is stuck around the "instantaneous" min- 
ima characterized by the condition + QFR = 0, i.e.. 



2Q(j>„ 



f 2VipC \ 

\ Pm ) 



< 1 



(19) 



where we used the fact that Vi is of the order of the 
squared of the present Hubble parameter Hq so that 
the potential ([5]) is responsible for the accelerated ex- 



2FH = -(1 - 6Q2)i;^^2 ^F + HF-p.^- |prad , (7) P^^^'^™ ^^^^^^^ ^ 



-2Q</.„ 



1 



{1-6Q^)F U + 3H^+^A +V^ + QFR = 0,{8) 



where R^^6{2H^ + H). 

In order to solve the background equations ([6])-([8|) nu- 
merically, we introduce the dimensionless variables 



xi 



We also define 




X3 



fluE= {l~6Q^)xl+xl + 2V6Qxi, (10) 

f^rad = xj , (11) 

n,n = 1 - (1 - 6Q^)xl -x\- 2V6Qxi - xl , (12) 



under the condition (jl9p . In this region the quan- 
tity |A| defined in Eq. (|14p is much larger than unity. 
The field value 10ml increases as the system enters the 
epoch of an accelerated expansion, which leads to the 
decrease of |A|. The matter-dominated epoch is real- 
ized by the instantaneous fixed point characterized by 
{xi,X2,X3) = (y6/(2A),[(3+2QA-6g2)/2A2]i/2^o) with 
rim = l-(3-12g2+7QA)/A2 ~ 1 and w^tf = -2Q/A ~ 
(because |A| ^ 1 in this regime). In the presence of 
the coupling Q there exists a de-Sitter point character- 
ized by (a;i,a:;2,a;3) = (0,1,0), = and WaS = -1, 
which corresponds to A = AQ. This solution is stable for 
dA/d0 < [23 and hence can be used for the late-time 
accelerated expansion. See Ref. [l^] for detailed analysis 
about the background cosmological evolution. 
The mass squared, Af^ = V^,^,, is given by 



which satisfy the relation f^DE + f^iad 
Eqs. and (O we find 



1. Using 



WiCpQ^il^pe 



-2Q<j, 



)(1 



e-2Q^)P-2e-2W.(20) 



Plugging the field value into Eq. ([20|) . we find 



H 



where A 



1-6Q 



2 r 



3 + Sx? - ixl + xl- GQ^xl 



hp 



l-p 



{2PpC)^/(^-p) 



Pm 



Vi 



(21) 



-2^/%Qxi 



3Q(Ax-2-4Q),(13) 



-V,^/V . For the potential ([5]) we have 

_ 2CpQe-2Q^(l_e-2Q'^')P-i 
- 1 - C(l - e-2Q*)P ■ 



(14) 



Since the energy density pm is much larger than Vi 
during the radiation and matter eras, we have that 
» Vi ~ i?o • The mass squared decreases to the 
order of Vi after the system enters the accelerated epoch. 
This evolution of the field mass leads to an interesting ob- 
servational signature in weak Icnsing observations, as we 
will see in subsequent sections. 
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III. WEAK LENSING 

Let us consider a perturbed metric about the flat 
FLRW background in the longitudinal gauge: 



d.s^ 



(1 + 2$)dt2 + a^{t){l - 2-^)Sijdx'dx^ , (22) 



where scalar metric perturbations $ and 4* do not co- 
incide with each other in the absence of an anisotropic 
stress. Matter density perturbations (S„i in the pressure- 
less matter contribute to the source term for the gravi- 
tational potentials $ and ^. The equation of 6m for the 
action ^ was derived in Ref . 12211 u nder an approxima- 
tion on sub- horizon scales [1, l35l. |37|| . Provided that the 
oscillating mode of the field perturbation S(j) does not 
dominate over the matter-induced mode at the initial 
stage of the matter era, we obtain the following approx- 
imate equation 

'in + 2H5,n - 47rGoffp™(5™ ~ , (23) 
where the effective gravitational "constant" is given by 
1 (fcVa2)(l + 2Q2)f + m2 



GeS 



SttF 



{k^/a'^)F + M^ 



(24) 



Here fc is a comoving wavenumber and is given in 
Eq. ([20|) for the potential ([5|). Using the derivative with 
respect to A^. Eq. (|23p can be written as 



diV2 



1 3 \ d<5,„ 
,2 - 2^-=^ diV 
(fcVa2)(l + 2Q2)i^ + M2 



2"""' (fc2/a2)i^ + il/2 
The gravitational potentials $ and ^E* satisfy 



Pm (fc2/a2)(l + 2Q2)F + M2 
2F (fc2/a2)F + A/2 ' 
p,„(fc2/a2)(i-2Q2)f + Af2 
'2F {k^/a^)F + M^ ' 



0. (25) 

5rn, (26) 

(27) 



In order to confront our model with weak lensing obser- 
vations, we define the so-called deflecting potential [s^ 



$wl = + * , 

together with the effective density field 

a 



Ooff 



(28) 



(29) 



where the subscript "0" represents the present values and 
we set ao = 1. Using the relation 



Prn — iFoH^VLmfl/ a , 

together with Eqs. ((261) and (H?]), we get 

^0 c 



3> 



wl 



Pm 



<5eff 



(30) 



(31) 



We write the angular position of a source to be 9s and 
the direction of weak lensing observation to be Oj. The 
deformation of the shape of galaxies is characterized by 
the amplification matrix A = d9s/d9j. The components 
of A are given by [H, [33 | 



A^i/ 



x'ix-x') 



X 



d^.-^^iix'O^x'W , (32) 



where x is the comoving radial distance satisfying the 
relation dx = —dt/a{t) along the geodesic. In terms of 
the redshift defined byz=l/a— l.we have that 



Xiz) 



H{z') 



(33) 



The convergence k and the shear 7 = (71,72) can be 
derived from the components of the 2x2 matrix A, as 

k^I-^TtA, {[A22-Aii]/2,Ai2) . (34) 

If we consider a redshift distribution p(x)dx of the source, 
the convergence is given by k{9) = / p(x)k(6', x)'^X- Us- 
ing Eqs. (Ull), ^ and ^ we obtain 



k{9) = ^H^nmfl I 9{X)X 



dx , (35) 



where xh is the maximum distance to the source and 

(36) 



g{x)^ p(x')^^dx'. 

Since the convergence is a function on the 2-sphere it 
can be expanded in the form k{9) = J k{i)e^^'^^, where 
£ = (^1,^2) with ii and £2 integers. Defining the power 
spectrum of the shear to be {k{£)k*{e')) ^ P„(f)5(2)(f- 
£'), one can show that the convergence has a same power 
spectrum as [2^ . It is given by Q 



Xh 



gix) 



a(x) 



£ ■ 
-iX 

X . 



dx- (37) 



In our scalar-tensor theory we have Pg^^^ ~ 
{Fo/F)'^Ps^ from Eq. ((3T|l . where Ps,„ is the matter 
power spectrum. In the following we assume that the 
sources are located at the distance Xs (corresponding to 
the redshift z^), which then gives p{x) = ^{x ~ Xs) and 
ff(x) = iXs ~ X)IXs- This leads to the following conver- 
gence spectrum 



Xs-xPq 



Xsa 



F 



Ps„ 



,X 



dx- 



(38) 

Let us consider the action ([2]) with the potential ([5]). 
In the deep matter era where the Ricci scalar R is much 
larger than Hq, we have M'^ / F k'^/a'^ and F ~ 1 for 
the wavenumber k relevant to the matter power spectrum 
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[22j . Since GeS ~ G in this regime from Eq. ([24]), the 
perturbations evolve in a standard way: Sm t^^"^ and 
<f>wi — constant. Meanwhile, at the late epoch of the 
matter era, the system can enter a stage characterized 
by the condition M^/F < P/a'^. Since GcH ~ (1 + 
2Q'^)/8ttF dm-ing this stage, the perturbations evolve in 
a non-standard way: 



The critical redshift Zk at M'^/F = k'^/a'^ can be esti- 
mated as 



Zk 



Hi 02(1 



2PpC 



- 1. 



(40) 

As long as Zk > 1 it is expected that the effect of modified 
gravity manifests itself in weak lensing observations. 

Since the evolution of perturbations is similar to that 
in the ACDM model at an early epoch characterized by 
the condition z ^ Zk, the deflecting potential <i>wi at late 
times is given by [i^ 

$wi(fc, a) = ^<i>^i{k, a,)T{k)^^^^ , (41) 
10 a 

where $wi(fc,ai) — 2$(fc,ai) corresponds to the initial 
deflecting potential generated during inflation, T(fc) is 
a transfer function that describes the epochs of horizon 
crossing and radiation/matter transition (50 ^ z < 10^), 
and D{k, a) is the growth function at late times defined 
by D{k,a)/a = <&wi(a)/'i'wi(a/) («/ corresponds to the 
scale factor at a redshift 1 <C z/ < 50). 

Since we are interested in the case where the transition 
redshift z^. is smaller than 50, we can use the standard 
transfer function of Bardeen et al. \4M ■ 



T{x) 



ln(l + 0.171a-) 
0.171a; 

+{ldM9xf + (0.490a;)4 



1.0 + 0.284a; + (1.18a;) 

-0.25 



(42) 



where x = k/k-EQ and /ceq = 0.073 r2,„_o/^^ Mpc ^. 

In the ACDM model the growth function during the 
matter-dominated epoch (f2m = 1) is scale-independent: 
D{k,a) = a} In our scalar-tensor model the mass 
squared given in Eq. pTjl evolves as cx 

^-2(2-p)/(i-p)^ .^yjjjch implies that the transition time 

at M^/F = k'^ /a? has a scale-dependence tk oc k 

[2^ . This leads to the scale-dependent growth of metric 

perturbations. 



^ Note that in the late-time accelerated epoch the growth of matter 
pcrtubations is no longer described by D(a) = a. 



Using Eqs. p9|) and ([4T|) we obtain the matter pertur- 
bation 5m at the redshift z < z/: 



5m{k, a) 



3 F 



10 Fq n,n,oHQ 



$wi(fc,a«)T(fc)i^(fc,a). (43) 



The initial power spectrum generated during inflation is 
P*^^ = 4|$|2 = {200TT^/9k^){k/Ho)'''-^Sjj, where n, is 
the spectral index and Sj^ is the amplitude of . Then 
the power spectrum, Pg^^ = |(5„ip, is given by 



P5,„ (fc, a) = 27r2 



F 



S%T^{k)D'^{k,a) 



From Eqs. and (gH) we get 



(44) 





X 



where 



^"EQ X 



dz , (45) 



(46) 



From Eq. ([55)1 the quantity X satisfies the differential 
equation AX/Az = 1/E{z). In the following we use the 
value Zs = 1 in our numerical simulations. 



IV. OBSERVATIONAL SIGNATURES OF 
MODIFIED GRAVITY 

When Q 7^ the evolution of 5m during the time- 
interval tk < t < t\ (where is the time at a = 0) 
is given by Eq. whereas Sm oc t^/s tj^c ACDM 

model (0 = 0). Hence, at time ^a, the power spectrum 
for Q exhibits a difference compared to the ACDM 
model 1231: 



pACDM(^^) 

where 



An(tA) 



\/25 + 48Q'^ -1 2 



oc fc'^"^*^) _ (47) 



(l-p)(v/25 + 48g2_5) 



P 



(48) 



In order to derive the difference An{to) at the present 
epoch, we need to solve perturbation equations numer- 
ically by the time to- However, as long as Zk is larger 
than the order of unity, the growth rate of Sm during the 
time-interval < t < to hardly depends on k for fixed 
Q. Hence it is expected that the analytic estimation ([48]) 
does not differ much from An(to) provided zj, ^ 1. 

We start integrating the background equations ([TB]) - 
P8[) from the deep matter era and identify the present 
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(a) Q=0 7, p=0.6 
(b) Q=-(1/6r^, p=0.6 

(c) ACDM, linear 
(d) ACDM, nonlinear 



o 




10" 



0.01 0.1 1 

k [h Mpc ■'] 

Figure 1: The matter power spectra Psmi^) the present 
epoch for (a) Q = 0.7, p = 0.6, C = 0.9, (b) Q = -1/V6, 
p = 0.6, C = 0.9, (c) the ACDM model, and (d) the ACDM 
model with a nonlinear halo-fitting (ag = 0.78 and shape 
parameter F — 0.2). The model parameters are Qm.o ~ 0.28, 

In 



= f and Sfj = 3.2 x fO" 



Ho = 3.34 X fO^" /iMpc' 
the cases (a) and (b) we start integrating Eqs. (|f 6|l - (|f 81) with 
initial conditions {xi,X2,X3) = (0, [(3+2QA-6Q^)/2A^]^/^, 0) 
and F - f = -fO"*. 



epoch by the condition = 0.28. Wc then run the 
code again from z = zj {< 50) to z = in order to solve 
the perturbation equations (^5)) and ([?T|l . Since we are 
considering the case in which Zk is smaller than z/, the 
initial conditions for matter perturbations are chosen to 
be ^ = S,n (i.e., those for the ACDM model). 

In Fig. [1] we plot the matter power spectra at the 
present epoch for (a) Q ^ 0.7, p — 0.6, C = 0.9, (b) 
Q = -1/V6, p ^ 0.6, C = 0.9, (c) the ACDM model, 
and (d) the ACDM model with a nonlinear halo-fitting 
H^. Since we do not take into account nonlinear effects 
in the cases (a)-(c), these results are trustable in the lin- 
ear regime k < 0.2/iMpc^^. 

In the case (b), which corresponds to f{R) gravity with 
n = 0.75 in the model the spectrum shows a devi- 
ation from the ACDM model for k > O.Ol/iMpc'^ On 
the scales k = O.Ol/iMpc"^ and k = O.l/iMpc'^ the crit- 
ical rcdshifts at A'P/F — k'^/a'^ are given by z^ ~ 2.995 
and Zk = 5.868, respectively. Numerically we find 
An{tn) = 0.017 and An(to) = 0.119 for k = O.Ol/iMpc"! 
and k ~ 0.1/iMpc~^ respectively, whereas the estimation 
(|48ll gives the value IS.n(tt\) = 0.088. Since Zfc decreases 
for smaller fc, the analytic estimation ([39|) obtained by 
using the condition z^ ^ 1 tends to be invalid on larger 
scales. This is the main reason of the discrepancy be- 
tween An(to) and An(tA) found for k < 0.1/iMpc~^. 
We checked that An{ta) approaches the analytic value 
An(tA) = 0.088 on smaller scales, e.g., An(io) = 0.089 
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(b) p=0.7, C=0.9 

(c) ACDM 
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Figure 2: The convergence power spectrum Pk{£) in f{R) 
gravity {Q — — 1/\/6) for the cases: (a) p = 0.5, C = 0.9 
and (b) p = 0.7, C = 0.9. We also show the spectrum in the 
ACDM model. Other model parameters are chosen similarly 
as in the case of Fig. [1] 



for k = 4.3/1 Mpc-i. 

For larger \Q\ the growth rate of dm increases in the 
regime z\ < z < Zk, which alters the shape of the matter 
power spectrum. In the case (a) of Fig. [T] we numerically 
find that An{to) = 0.323 on the scale k = 0.1/iMpc"\ 
while the estimation (|48|) gives An(<A) = 0.231. Again 
this analytic estimation is in a better agreement with 
An{to) on smaller scales, e.g., An(to) = 0.244 for k = 
4.3/iMpc~^. In Fig. [T] we also show the matter power 
spectrum in the ACDM model derived by using the non- 
linear halo-fit [i^. This gives rise to an enhancement 
of the power in the nonlinear regime (fc > 0.2/iMpc~^). 
The spectrum in the case (a) exhibits a significant differ- 
ence compared to this halo-fit ACDM spectrum even for 
k < 0.2/iMpc~^, which implies that our linear analysis is 
enough to place stringent constraints on model parame- 
ters Q and p from observations of galaxy clustering. 

Let us next proceed to the convergence power spectrum 
of weak lensing. Compared to the matter power spec- 
trum the wavenumber k is replaced by fc = (/x- the 
deep matter era the evolution of the Hubble parameter 
can be approximated as H'^{z) ~ HQilmfii^ + z)^, which 
gives X — 2/(-ffo^^,no) ^ constant. Hence the time ti at 

AP/F = (^/x)^/a^ has an ^-dependence ti oc £~ , 
provided this transition occurs at the redshift ^ 1. 

Since $wi — constant for tj < t < ti and $wi oc 
^(v/25+48Q^-5)/6 f^j. tg < t < tA, wc have that 



^wi(za) 
$wi(z/) 



tt J 



(V25+48Q2_5)/6 



(49) 
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(a) Q=1, C=0.9 
(b) Q=0.5, C=0.9 
(c) ACDM 
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Figure 3: The convergence power spectrum Pk{£) for p = 0.7 
with two cases: (a) g = 1, C = 0.9 and (b) Q = 0.5, C = 0.9 
together with the ACDM spectrum. Other model parameters 
are cliosen similarly as in the case of Fig. [T] 



find that ze = 1.386 for £ = 5. In this case the system 
enters tlie non-standard regime (z < zt) before entering 
tlie epoch of an accelerated expansion (z < z\ ^ 1), 
which leads to the amplification of <i>wi. This changes the 
total amplitude of Pk{£) relative to the ACDM model. 
The numerical value of An at £ — 200 is found to be 
An = 0.084 for p = 0.5. Since An increases for smaller 
p, this information is useful to place a lower bound on p 
in f{R) gravity from weak lensing observations. 

In Fig. [3] the convergence spectrum for p = 0.7 is plot- 
ted for two different values of Q together with the ACDM 
spectrum. We note that the transition redshift Zi de- 
creases for larger \Q\, see Eq. PD|) . Hence the deviation 
from the ACDM model is insignificant for small £, unless 
we choose smaller values of p. However the spectrum is 
strongly modified for £ > 10 with the increase of |Q|. 
The numerical values of An at £ = 200 are found to be 
An = 0.084 and An = 0.311 for Q = 0.5 and (5=1, 
respectively. Hence it should be possible to derive an up- 
per bound on the strength of the coupling Q by using 
observational data of weak lensing. 



CONCLUSIONS 



As long as z^ 1, the evolution of $wi during the time- 
interval tA < t < to is almost independent of £ for a fixed 
value of Q. Then we obtain the following .^-dependence 
for < z < ZA ^ z^: 



^wl(-g) 
$wl(z/) 



cx £- 



(l-p)(\/25 + 48Q^- 



(50) 



From Eq. p5)) this leads to a difference of the spectral in- 
dex of the convergence spectrum compared to the ACDM 
model: 



PACDM(£) 



(51) 



where An is the same as An(tA) given in Eq. We 
caution again that the estimation (j5ip is valid for z^ 3> 1 . 

In Fig. [2] we plot the convergence spectrum in f{R) 
gravity for two different values of p together with the 
ACDM spectrum. We focus on the linear regime charac- 
terized hy £ < 200. Since the ACDM model corresponds 
to the limit n — > oo in Eq. ([3]), the power p = 2n/ (2n+ 1) 
approaches 1 in this limit. The deviation from the ACDM 
model becomes important for smaller p away from 1 . 

When p = 0.7, for example. Fig. [2] shows that such a 
deviation becomes significant for £ ^ 10. Numerically we 
get An = 0.056 at £ ~ 200, which is slightly smaller than 
the analytic value An = 0.068 estimated by Eq. ([5T|) . 
The main reason for this difference is that the critical 
redshift z^ = 3.258 at £ = 200 is not very much larger 
than unity. 

When p = 0.5 the deflecting potential <f>wi is amplified 
even for small £ {< 10) , which is associated with the fact 
that zi is greater than 1 even ior £ > 2. For example we 



We have discussed the signature of modified gravity 
in weak lensing observations. Our model is described 
by the action ^ with a constant coupling Q, which is 
equivalent to Brans-Dicke theory with a field potential 
V. This theory includes f{R) gravity as a special case 
(Q = -l/\/6). The scalar-field potential V{cj)) can be 
designed to satisfy local gravity constraints through a 
chameleon mechanism. The representative potential that 
satisfies LGC is given in Eq. ([5]), which is motivated by 
viable ,f{R) models proposed by Hu and Sawicki [l^ and 
by Starobinsky |1Q] . Note that most of past works in 
scalar-tensor dark energy models restricted the analysis 
in the small coupling region {\Q\ < lO^'^). In this paper 
we focused on the large \Q\ region in which a signifi- 
cant difference from the ACDM model can be expected 
in weak lensing observations. 

Cosniologically these models can show deviations from 
the ACDM model at late epochs of the matter-dominated 
era. The growth rate of matter density perturbations gets 
larger for redshifts smaller than a critical value z^. Since 
Zfc increases for larger k, the matter power spectrum is 
subject to change on smaller scales. We evaluated the 
matter power spectrum Pg^ (fc) numerically and showed 
that the spectral index and the amplitude of Pg^^ (fc) can 
be significantly modified for larger values of \Q\. 

The non-standard evolution of matter perturbations 
affects the convergence power spectrum Pk{£) of weak 
lensing. As long as the transition redshift zi is larger 
than the order of unity, one can estimate the difference 
An of spectral indices between modified gravity and the 
ACDM cosmology to be An ~ (1 - p)(^25 + 48Q2 _ 
5)/(4 — p) with < p < 1. In f{R) gravity the parameter 
n for the model ^ is linked with the parameter p via 
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the relation p = 2n/{2n + 1). The hmit p — > 1 (i.e., 
n ^ oo) corresponds to the ACDM model, in which case 
we have An 0. The difference of the convergence 
spectrum relative to the ACDM case is significant for p 
away from 1. As seen in Fig. [2] (which corresponds to the 
case Q = — 1/ \/6) , the spectral index and the amplitude 
of -Pk(^) are modified for smaller values of p. 

If we take larger values of \Q\, the convergence spec- 
trum deviates from that in the ACDM model more sig- 
nificantly. This situation is clearly seen in the numerical 
simulation of Fig. [3] It should be possible to place strong 
observational constraints on the parameters Q and p by 



using observational data of weak lensing and the mat- 
ter power spectrum, which we leave for future work. We 
hope that some signatures of modified gravity can be de- 
tected in future high-precision observations to reveal the 
origin of dark energy. 
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